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ABSTRACT 


A single linear vertical passive array is used in the 
'SOFAR' channel to determine the depth of a single under- 
water source at a constant range. The phas¢ and amplitudes 
weights applied to the array are determined by the linear 
minimum variance estimation technique. The resulting beam 


pattern is compared to the conventional time demain bean- 


ct 


Gobo 


former. It was found that the linear minimum variances 
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mation «¢chnigue of amplitude shading and phase weighting 


was significantly superior to the conventional beamformer. 
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I. INTRODUCTION 

There are many sclutions to the problem of voredicting 
the location of an underwater energy scurce. Oné conmon 
soluzicn is the use cf & passive hydrophone which detects 
the pressure waves radiating from the source. The hydro- 
phone sensor is assumed to be omnidirectional and therefore 
incapable c= estimating direction. io provide direccicn- 
ality a series cf sensors are placed in a row to fcrma 


passive linear array. 


7) 


A familiar methed of determining directionality i 
time-domain beamforming. In this principle, it is assumed 
*hat the source is far enough away so that the pressure wave 
appears to be a plane wave when viewed at the site of the 
receiving array. : 

mais 2 Set Of: time delays are calculated for any direc- 
tion cf signal arrival, which, when applied to the receiv 
outputs causes them to be in phase and to reinforcs when 
summed. -The resultant angular response to signals arriving 
from other than the nominated direction is «hen a function 
ct the array geometry, relative to the signal wavelength, 
and any weighting factors which have been applied to «he 
recéiver outputs. The effect is to generate a main 
receiving beam in the desired direction, with a saries of 
undesired subsidiary sidelobes whose magnitude can be 
controlled to some extent by EUewmcnO LCG Gf a SUicable array 
geometry and the use cf amplitude weightings on the rectiver 
Sucputs. 

In crder to determine location, three or more such 


arrays separated by a known amount may be used. 


Hi 
ct 
rw 


This study is ccncerned with a single linear ve 


c 
passive array and the determination of the depth of a single 





underwater source. The analysis is based on the fcllcwing 
assumpticns: 
e The underwater source is ¢mitting continuously and ata 


monochrcematic frequency. 


e Both the source and receiving hydrophones are stationary 
in space causing a constant range. 

e The range is sufficiently long so that the channsi is 
filled with R-R (refracted-refractsd) rays. 

aeese £5 7 distortion introduced in the prcepagating 


medium so that the signals received at each sensor are 


identical except for constant délays. 


e The source Signal and noise aré 
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staticnary gaussian tandom processes. 


Seine =pecd of sound profile is triangular and symmetric 
with the deep scund channel axis a= 1000 meters. The 
velocity gradient is -0.017 meters/meterysec above 100600 
meters and +0.017 méters/meter/sec below 1000 meters. 
This profile gives a speed of sound at the surface of 


1500 meters/sec. 


Bele screed cf sound profile is constan= in the herizcntal 


plane. 
e Only R-R rays are considered. All ocher rays have 
Suctsczent less that “heir eff]est is reqligible. 


As opposed to conventional <cim2-domairz 


a 
ao 
wo 
ip 


SOrn2ig, chilis 


nl 
e fronts at the 


study makes no assumption cf planar wav 

teceivezt site. Therefore the time delays applied to each 
rec2iver will not, in general, be a ilinsar functicn of 
depth. 


Since i+ is desired +o determine whether or not there is 


a scurce present at a specific depth the result will bea 
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binary decision. ADM weil andicate siqnal source 
present; a "0" will indicate signal source not present. For 
the constént range there will be "N" test depths investi- 
gated for the signal source. The number of hydrophon 

the vertical array will be "L". 

Fer a single source at a given depth, the travel time is 
calculated from the depth to each hydrophone. This travel 
time is ccenverted intc @ phase delay for each hydrophone so 
that after summation from all hydrophones a maximum output 
is achieved. This cutput is then passed thrcugh a squaring 
@=v-ce, an integrator, and a threshold and flip fler device 
comaga.ve a 1" binary output. If the signal source is ata 

ifferent depth and the same previous phase delays are used 

for each hydrophone then the output will be somawhat less 
than the previous maximum. The difference in depth required 
to achisve a"0O" binary output is the depth resolution of 
the system. 

The travel times for each hydrophone are calculated for 
each ci the "N" source depths =o be considered. VEN Sw as 
crdinarily be much greater than "L" so that the system will 
be overdetermined. An overdetarmined system is one in which 
there are more equaticns than unknowns. The objective then 
is to calculate the phase angle and amplitude weight for 
each hydrophene so that a determination can be made indi- 
cating the presence or absence of a signal scurce at a given 
depth. 

The method used to calculate the phase and amplitude 
weights is the linear minimum variance ¢stimation technique. 
Iinear minimum variance estimators are optimum when compared 
with all cther estimators for gaussian problems. The method 
is directly applicable to overdetermin2ed systems. 

The output of the summer is calculated using the linsar 
Minimum variance amplitude and phase angles assuming a 


source at one of the "N" source depths and no source at the 


an 





ethers. MTh+ calculation is repeated for each of the dertns. 
The result, when pictted against source depth, wiili be 
oas the “beam pattern" of the array 
(Although similar, it shouid not be interpreted 2 
the angular response of an array aS in the conventional 
definiticn cf a beam pattern. The conventional definition 
loses much of its utility when the wavefronts are not 
Eeoder. ) Ideally *he beam pattern will be maximum at the 
desired depth and very small at ail other depths so that the 
binary "1" decision will be made for a source at the desired 
depth, anda "0" decision for sources az all others. This 
beam pattern is compared with the depth bean pattern of the 
conventional time-domain beamformer mentioned above. The 
beepcee cf this *hesis is te determine, as an initial inves- 
tigation, whether the linear minimum variarc? estimation 
technique, when applied to a linear vertical array, is 
useful in depth discrimination at long rarges ir 2 'SOFAR' 
type sound channel. 
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A. RAY ACOUSTICS 


The propagation cf sound in an elastic medium can be 
described mathematically by solutions of the wave equation 
using the appropriate boundary and medium conditions fora 
particular problem. The wave equation relating the accustic 
Pees=ure 'p' to the coordinates 'x','ty','z', and the time 


'c', may ke written as 


doe 


He? de 


2 Be ne 
dp eG dsp 


= + 
7 





AO fe (2.1) 
whezte 'c' is a quantity that has the general significance of 
sound velocity and may vary with the coordinates. 

One may approximate the solution of the wave equation 
using .ray theory: its bedy of results and conclusions is 
called ray acoustics. 

Officer [{Ref. 1] describes the ray solu+ion as a 
complet? soluticn +c any particular propagation opreblem 
Beeoin the validity cf the approximation of the Eikonal 
equation to the wave equation. For these approximations to 
ke valid reither the amplitude of the wave nor the speed of 
scund can change aprreciably in distances compatabls to a 
wavelength. 

Thus the path of a ray through a medium in which the 
Seeeopo. Sound varies with depth can be caiculated by the 


ow 


pplicaticn of Snell's law 


cos@/c = Ve =a constant for any one ray (252) 
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where '€' ais the angie cr depression made with hs 
herizontal at a depth where the speed of sound is 'c', and 
"Co ' is th2 speed at a depth (real or extrapolazed) 

the ray would become herizontal. 

In ame2dium in which the velocity of sound changes 
linearly with depth +he sound rays can be shown to be arcs 
Semcarcies, that is, to have a constant radius of curvature. 
Kinsler ez al. [Ref. 2] give a simpls and heuristic demons- 
Beeson Of the circularity of rays ina medium with a lLinsar 
sound speed gradient ‘'g'. Miescentsr Of the cGirele which 
creates the arc lies at a depth where th2 sound speed extra- 
polates to zero. To understand this, considar a porticn of 
Pueeoye Path with a local radius of curvature 'R', as illus- 
feeated in Figure 2.1. Since the gradient 'g' for this case 


is 
g=Ac/Az = (C,-C,) / (d4-d4) = (Cy -C) /R (C988) -COS8 >) (2,3) 
where 'ac' is the change in sound speed and ‘,jz' is the 


Change in depth. It can be seen that the radius of curva- 


ture is given by 


R = ae. /9 = =¢c/(9 Ccss) (2.4) 
The ray path is therefore a circle when 'g! is constant 
because 'R' is then conszant. The Cenc®er of Curvature cf a 


circle lies at the depth where @ is 90 degrees, which corre- 
sponds to c=0. Hor the Situa.2on in Figure 2.1 the spsed 
gradient is negative so that 'R' is positive. If the speed 
gradient were positive 'R' would be negative, and the path 
would curve upward. 

Once the radius cf curvature of éach segment cf a path 
is known the actual path can be traced graphically or 
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Velocity 
protile 





Figure 2.1 ircular Bay Pati 


1: 





computed. If the initial angle of depression of a ray is 
61, then by referring tc the geometry of Figure 2.1 the 
changés in beth range and depth are 


Ar= c, (sin (91) -sin (@2))/(g cos (91) ) (225) 


AZ= Cz (cos (82) -cos (61) )/(g cos (81) ) (2.6) 


The sign convention for these equations is: downward, «=c che 
ma ght , and depressicn angles below <=he horizcntai axis 
positiv:. 

The symmetric triangular sound spted profil2 assumed in 
the introduction is similar to spezd profiles encountered in 
the deep sound channel, sometimes called the SOFAR channel. 
The velocity minimum which occurs at the axis of zhe sound 
channel causes the sea to act like a kind of lens; abeve and 
below the minimum, the velocity gradient continually bends 
the scund rays soward the depth of minimum velocity. A 
pertion of the power radiated by a source in the deep sound 
channel acccrdingly remains within the channel and encoun- 
ters ne accustic losses by reflection from the surface and 
bottcm. These rays are called R-R (refracted-refracted) and 
Since tkey have very low transmission loss, very long rangés 
can be obtained from a source of moderate acoustic power 
Sie put. Thus an energy source at a specific depth will 
propagate energy in all directions but cnly che direction 
which is tceward the receiving array and for which the rays 
are R-R is of interest. To determine the range of derres- 
sion angles which will yield R-R rays Snell's law is used to 


give 
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@max = arc cos (cl/ (ci-(g d))) (oa) 


where ‘O6max' is in radians and 'cl' is the speed of sound at 
foesscurce depth 'd'. 

An @xampl= is given in Figure 2.2 of a single ray trace 
propagatirg in the SOFAR channel to show how aquaticns 2.5 
and 2.6 are used in determining range and depth. The ray 
frets Exoken up into arcs of citcles as shown in the 
maguce, and then by paying close attention to the previously 
defined sign conventicns, the chang? in range and depth is 
mound. Be nguMone spec in2¢, for arc 1, @1 and @2 are both 
Pereecoy=:- £05 arc 2, 61 25 positives and @2 is zero; for are 
3, Simeec Ze>0 and 62 is negative: and for arc 4, #61 is 
negative and @2 is zero. By keeping @ running total of all 
the déepcth and range changes it is possible ts determine che 
Seeaenorazon-al distanc2 travelisd and the depth at that 
distance. 

For the speed of scund profile assumed in the inztroduc- 
maomee ) computer generated tay plot is shown in Figure 2.3 
for a source depth ci 300 meters. AS each ray propagax<es 
out from the source the triangular channel becomes filled 
Wath sound. -Il£ a receiving hydrophone is placed a gréa 
distance away, a number cf refractsd propagation paths will 
exist, each having a differen: travel tine and ctossing <zhe 
Omens! axis at different intervals. The path with the 
greatest excursion from the axis will have the shortest 
travel «ime. 

Officer [Ref. 1] shows that the travel ctime 't' of a 


Maeeeewn ch 25 am arc of a circle, is given by 


82 ae 
S cosé (2.9) 
8 





peas che travel tfime for each arc is 
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8 iors tan((17/4)+(62/2)) (229) 
ne tan((7/4)+(8;/2 


Equation 2.9, applied using the same convention as equations 
maoeartd 2.6, determines the total travel time of a ray in 
the deep sound channel. 


Be. ARRAY MODEL 


The receiving linear vertical array is presumed <o 
consist of 'L' hydrcphones aS shown in Figure 2.4. cis 
assumed that the source iS emitting energy at a constant 
foequency, ‘ft, and amplitudes, ‘A’, regardless of the depth. 
The source signal at the source is Aexp(j2™ fr). The 


inherent received signal at the first hvydropnonse is 


x1(t) = Atexp (j2TE (t-t1)) (22710) 
where 't1' is the travel time from the energy source to the 
first hydrophone and ‘'At' is the amplitude of the signal at 
the range cf the array. After passage «hrough the amplitude 


weight ‘ai’ and a pkase delay of ‘'T1', the signal on che 


first hydrophone at the input to the summex is 


yitt)=A'alxi(t-11)=A'alexp(j27f(t-t1l-t1)) (2.11) 


A time delay, for monechrematic siqnais corresponds =o a 


phase shift 


@ = 2nft . Does 


where '@' is the phase shift in radians and '*' is the time 


delay in seconds. Thus éguation 2.11 can be written as 
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y1(t) =Atatexp(4(2 £t-$1-@1)) (es) 


where ')1=20ic1! is the phase delay due to the travel time 

Meemetie source tc «he first hydrophene and '6@1' is the 

phase delay in the receiver on the firs hydrophone. 
Combining all the hydrepkhones in the arzay in a sumzer 


gives 25 an expressicn for the array output 


L 
y= Bak eRe (LATE fy 8) ) (2.14) 
where oy! represents the phase dslay due zo the travel time 
Peemecais source to the "“"k th" hydrophone and "O,' is the 


3 


phase delay in the receiver on the "k th" hydrophone. ie 
the amplitude of the energy source is normalized by seczting 
At=1, and equation 2.14 is written in ‘+erms of real and 


imaginary ccmponents, we have 


ay. cos(-$,-6, )+jsin(-$,-8,) expo j ener) (2. 15) 


When an energy source is at tha "q th" depth, we wish to 
hav2 each receiving hydrephcne's phase delay cancel out “che 
effect of the travel time Irom th2 source to it, such that 
tab EK! is equal to a multipls of '27', Tiss wit pas al) 
Signals into the summer in phase and thus maximiz? the 
Signal gain for a source at the "q th" depth. From equation 
eeilS ¢ 


oe © = L (source present at q) (2. 16) 


e (2. 17) 
) (aq SEB On eee = 0 (source present at q) ° 
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Note that the first subscripz on z=he phése angle indicates 
the recéivizrg hydrophcne and the second subscript indicates 
the depth of the sourcs. Thus "Dieg 'weulda  2pdicaze the 
PHase shift relating to the travel time from the "gq th" test 
@epen to che "k th"' hydrophone in the receiving array. 

It is desirable fer '‘'Y(t)' to b32 a minimum value for 
sources at cther than the depth being investigated. Thus 
Meteedeh Of tne other 'N=-1' depths 'Y¥(t)' is set to Zero. 
This gives 'N-1' equations for the real terms of "Y¥ (+)! set 


ele) AS wae) 


" : 
j ag Cos (~bm-Oxq) =0 (source absent; mAq) (2.18) 

kel 

and 'N-1' equations for the imaginary cerms of 'Y(t)' set to 


ZEro 


21g S20 Pm” Ox) 79 (scurce absent; m#aq) (2. 19) 


By using elementary trigoncmetric identities, equation 
2.16 (real terms with source present at "q <“«h" depth) 
kecomes 

L 
akg [eos xq) EOS (8kq)-Sin (Pkg) sia (2xq)] =L (2. 20) 


Equation 2.18 (real terms with scurc2 absent for each of the 


ther 'N-1' depths) becomes 


ue je 


r Big [228 hn? 208 Pg? “FERC Fgg) NCO) = m=1,2,...N; m#G (2.21) 


Bamgeacton 2.17 (imaginary terms with scurce present at the "q 
th" depth) becomes 
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L 
-a (sin cos (6 +ccs sin (@ =0 Dee 
Ss ve Ui On ee 2,0) (2. 22) 
and equation 2.19 (imaginary <erms with sourcs absent for 


each of the other 'N-1' depths) beccm:3s 


L 
sa Ma [P2497 €05(83.4) 4008 (4,,)5iN(8, SOmmim= 2) oN. mia (2833) 


Thus, Pemo are a tocal of "2N" equations with '2L' 

unkncewns. 
In ordez> to simplify, w2 put these '2N' equati 
*rarily the real terms a 


re) 
Maer ox. fCrom. Arbi (we WGle = qe) jeniere 
n N 


'N' equations and tke imaginary t2rms the second 'N' esqua- 
tions. The first real and first imaginary equation is at 
the lewest (shallowest) source depth and equations increase 
in order after that until the last real and last imaginary 
equaticn correspond tc a source at the deepest depth. The 


resultant matrix equation becomes: 


0 coséyy cosh y rere cos$) 4 
O{- cosé), cost. aious cosé, 5 


-sino,, -sino,, 











-sin,, ~sino,, 









0 

L iG SG 660 Bes eon Sea aG wets 
Q 

0 COS ba 4 COSbyy ee cos, y 

0 -sino,, -sino,, --sing, 4 -cosos, -cos$,, 

0 -sind,, -sind,, --sing, » -cos,, -cos$,, 

0 -sindsy -sing,, --sind 


-cos¢ -cos¢ 


Simplifying further, equation 2.24 becomes: 


ea 








(225) 





were a ', eho vappro pr. sce 
submatrices. 

mien, ty noting that the sultiplicaticn of each elenen= 
of a cclumn by the same nonzero constant doesn't affect the 


soluticn, equation 2.25 becomes: 





(2. 26) 
=f a S 





Mmiaetay tecrting Z, A, and § represent the matrices in ¢equa- 


tion 2.26, we obtain 


Z = he (2227) 


Seees-meee Mder2x iS denoted Ey 4a capitelized underlined 
letter. 

In summary, 2Z is the desired response of the vertical 
azray to the 'N' source test depths. A represents known 
travel times from each of the 'N' scurce depths to each of 
the 'L' receiving hydrophones. 2 is unknown. It is the 
phase and amplitude weighting which must be applied *c the 
vertical Seayerneecrde= ~<G tealize 2. for come d = 0S 02 21" 


unknowns. 
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BaWatacre 2.27 represents '2N' equations. Sines =his 
system of equations is overdetermined (N>L) an #xac 

mean, dees not exist. In order to mak? the best esczin ° 
@ for the desired response, the linear minimum variance 


estimaticn technique is used. 


C. LINEAR MINIMUM VARIANCE METHOD 


Equation 2.27 tepresents a noise frie environmen. ee 
noisé w2re present it would becomes 

cae BO +0 (2. 28) 
Beem 'n* a "2A" slement column vector. This represents the 


noise at each source depth. ‘'Z' is a linear function of ‘et 
@ad is called the oberservaticn. Wives rae Tonk ne dula- 
tion or cksérvation matrix which is known, and ‘e!' is the 
"2L" element random parameter vector which is to be esti- 
mated. Assume “<he first and second moments of '9' and 'n't 


are given by 


E(8) = Us, Var(®) = Vo (2.29) 
and 

E(r) = 0 Var(n) = Vn (2.30) 
where 'E' represents expectation or first moment and ‘Var? 
represents the covariance matrix. It is assumed that the 


parameter *'€' and the noise 'n’ are uncorrelated. 
A testriction imposed is that the estimate must be a 


weighted linear combination of the observations: 
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@ = b+tEZ ean 


where '*' indicates estimate. Pics Objective Le <tco select 
to and *B' in order te minimize the 2rror variance. Suen 
an estimator is the linear minimum variance estimator; it is 
the kest, in the sense cf minimum-error-variance linear 
estimators. 

Another restricticn is that the estimator be unbiesed; 
in otherwords it is required that ths expected value of the 
estimatcr oe is equal to the expactesd value of the paran- 


eter ‘et. Thus, 


B (8) =E+BE(z) = E(@) = y, (2732) 
yielding 
Pala -EAU, (2.33) 


Peetu ing this cesult in equaticn 2.31 gives for the 


unbiased linear estimator 


) (2. 34) 


Note that since the estimator is unbiased, the estimation 
error "978-6 ' is zero mean. Bile hede Se Pete c Selec: 
meh2n crder tO Minimize the error variancs. However, «his 
Cptimization problem is ill-defined because the error vari- 
ance is a matrix. Therefore in order to introduce a scalar 
goodness measure the sum of the variances of gach component 
peo’ is minimized. This is the sum of the main diagonal 
terms of the covariance matrix and is defined as the trace 


cf the matrix. 
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2N 
Sapo ear(e)) = ) Vart(a.)_) (ea) 
le, el ie) IA 


where Mee’ Syndicates trace. ‘pe eSecnen seleceed ste 


Mepamize the trace of ths errcr variance, cr 


min tz (Var(@_)) = min tr (E(@_82)) (2. 36) 
where "lt indicates the transpose. The following problen is 


then cbhtained by substituting equation 2.34 into equation 
Ze 30s 


. ae = = 7 = 
ae ip) Beer ee)) men tr CEC (9 uw, -3¢Z—Au, )(9-u,3¢(Z-Au,) )) (2. 37) 
momewell known [Ref. 3] that equation 2.37 is minimized 


when 


Cov (8,2)-BVar (Z) =0 ie 5) 
where 'Cev(@,Z)' is the ccvariancs matrix of the unknown 
parameters and the cbservaticns. Deter g Sate) Obein1n 
Balter by ip, te de=2t alse 


B=Cov (8,2) (Var(z)) ~+ (2. 39) 
aminimum is achieved for the sum of the squares of the 
errors. 


Meeng equattcn 2.28 for 'Z', PiouCOVmesLeowor No anid 
'Z' becomes 


Cov (9,2) =Cov (@,AG+n)=V,a" (2.40) 


orcad ie sabe Uncorrelated- Ths varlance of 'Z* is 
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Var (Z) =Vaz @@+n) =AV 


atey 2. 41 
0= =n hei 


Substituting these into equaticn 2.39 gives 


2 = VpA (AVA +¥,)7™ (2.42) 


and the linear minimum vériance estimator is 


x it q -1 

= = 4 
S8iuy = He tVgA CAVGA +V) ~(Z-Au,) Bon) 
By utilizing a matrix inversion lamma (Ref. 3] equation 


ze43 becemes 


A 


a jie 
Boy - ‘AY, 


TANG)“ T CAT yt z4v5*u,) (2. 44) 
The advantage of this equation over equation 2.43 i 
size cf the matrix tc be inverted. In equation 2.4 
Matrix has dimensionality ‘2N' whil in equation 2.4 
BPemenaonality is only '2L'. Thus the advantages of the 
linear variance ¢Stimator are the ease with which they are 
derived, «he mathematical tractability of che linear forn, 
and the mirimum amount of stochastic information required 
for develcpmen<. ANeaintereStingechareacreris=ic is that the 
linear minimum variance estimate is the erthegonal prcjec- 
meen cf "e' onto *he space spanned by the observation "Z'. 
Eecause of these factors this e¢stimatcer is 2 popuiar form 
for estimating unknowrs in cverdetermined equations. 

For this thesis it is assumed that the neise samples are 
uncorrelated and identically distribut2d so thas: 


2 
wo=-oLl (2.45) 
n 
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No previous knowledge is assumed abou~ '6@'. This implies an 


infinite variance matrix which is tepresented as: 


ey a and Ug = (2.46) 


The linear minimum variance ¢stimazta given by equation 2.44 


is then 


“Aw 


eri T 
= 247 
eyg(A AZ (2.47) 


By determinin Soy the phase and amplitude weights are 


found for a signal scurc2 on «he ‘'q th! depth. Recall that 





34 Sys SS Lg) 
30 Bq 9) 24 
3 2G = a, cos3 
SSeM Vi a — La L 
ret aoe 1g 
; (2.48) 
Bor a Goon a 


Upon solving, this equation gives for the phase delay 


nw “Aw aw 


= 2 =e m=1%* 5 
ae arctan (9, //@ ) whe m eo L (2.499) 


- 


and amplitude weigh 


a =8@ /cos(9@ ) where m=1 to L (2.50) 
mq m mq 


When these amplitude weights and phas2 delays are applied to 
*he vertical linear array a resulting beam pattarn is formed 


which in the absence of noiss is: 


7a 














2 =A (2c) 
The resulting beam pattern 'Z2 ' can then be compared with 
the desired beam pattern 'Z', as wall as with a conventional 


beam pattern 'Z'' obtained using linear phase shifts across 
the array apertur2 selected to "stsar" the arvtay to the 
dominanz arrival angle for the selacted source depth. 
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III. EXPERI MENTA 


PROCEDURE 


A. BASIC ASSUMPTIONS 


The speed of sound profile given in the Introduction was 
used to test the techinique, with 'L', the number cf source 
d2pths chosen as 20. The scurce depths are at 220 meters 
and every 20 meters thereafter +09 600 meters inclusive. 
Nt, the number of hydreopheres used in the vertical array 
was chosen as 5. The hydrophores are placed at depths of 
100, 200, 3200, 400, and 500 meters. Thus, theres are 40 
equations, two for each source depth, and 10 unknowns, ctwe 
for each hydrophcne. A range of 200 kilometers (km) assures 
that the deep scund channel is filled with sound over the 
aperture of the vertical linear array. A frequency of 100 
Hz provides good resclution in the beam pattern without 
introducing alias mainlobes at the selected range across the 


depths cf investigaticn. 


B. "A* HATRIX CALCULATION 


Since the gradient of the sound velocity, Wg “is 
constant in the area of the source depths, and the speed of 
sound at the surface is known, equation 2.3 is used to solve 
for 'c2', the spzed of sound a* the source depth. Equation 
2.7 is used for each source depth to calculate the maxinun 
initial depression argle which yields R-R rays. In deter- 
Mining sravel time, cnly depression angles from each scurce 
which are pesitive (dcwnward) and yield R-R cays are used. 

Beginning with the first source depth of 220 meters, an 
initial depressicn angle of 0.0 degrees is selected. The 
ray path is then calculated using aquatiors 2.5 and 2.6 and 


broken into a series cf arcs as in figure 2.2. Equation 2.9 


aD 





is used to determine the travel time for each arc in the 
same mMannér. Each are's herizontal range, travel tim2, and 
depth are sumned. When the summed horizontal rang¢ teaches 
200 km the summation process ends. The travel time and 
depth of the ray at this horizontal range is then knewn. 
The same procedure is repeated for an initial depression 
angle of 0.1 degrees and every increment of 0.1 déegrses 
thereafter until the maximunr depression from equation 2.7 is 
reached. The final ray path is at this maximum depression 
angle. 

The same procedure is repeatec for each of the other 19 
source depths. 

Thus, for each initial angl2 from each source d-pth 
there is a ray which has a travel time and a depth when it 
reaches tke horizontal range of 200 ka. Sin¢ee at is*.he 
profile of the sound préssure wave which impinges on the 
vertical array which is of importance, a constan= can be 
Subtraczed from these calculated travel times. Toss 
constant is selected to be the travel time for the source 
depth of 220 meters which has an initial depressicn angle of 
Q degrees. It is subtracted from each of the travel times 
making =the resultant travel times relative with respect to 
the ray which has a 0 degree depression angle from the 220 
meter depth. The pregram and its listing which calculates 
the relative travel times and the dépths of these rays at 
the horizontal range of 200 km is given in Appendix A. 

A plet cf the relative travel times versus depth for the 
220 meter source depth is shown in Figure 3.1. Figure 3.2 
displays the plot for the 380 meter source depth. Negative 
relative travel times in the plots indicate +hat the overall 
travel time is less than the reference. These rays arriv¢ 
at the 200 km horizontal distance before the reference ray. 

Since the receiving hydrophones are at set vertical 
positions (100, 200, 300, 400, and 500 meters), an 
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interpolation is done to determine relative travel times <o 
them frem ¢ach scurce depth. DNesantsnpolat=-on program and 
its listing is in Appendix 8. Sometimes more than one R-R 
ray travels from the source depth to a hydrophone. when 
this cccurs, the ray which arrives first is used in cuhe 
calculaticn of relative travel time to that hydrophone. 

Equation 2.12 determines the phase shift relating «o the 
relative travel times. The "At matrix is formed by taking 
the apprepriate sine and cosine values as in equation 2.24. 
wae “AY matrix is "40 by 10°. 


C. 2" MATRIX 


Bemeorring tO equation 2.26, the ‘Z' matrix is a ‘40 by 
1* column vector. It is the desired beam pattern. The 
bottom 20 rews give the imaginary terms and are set to zero. 
The top 20 rows represent the value of the real terms at 


each source depth. Therefore each of the top 20 rows is ser 


to zero except for the row containing the source. It is set 


Somle Fcr sxample, if the source is at 220 meters then only 
the tcp rew is set tcl. If th? source is at 380 nevers 


then only the nirth row is set to 1. 


D. RESULTING BEAM PATTERN 
1. Using the Linear 


"eM ' is calculated using equation 2.47. The 
festering team pattern 'Z2* is calculated using ¢quazion 
Zee | « The prcgram which calculates the 'At matrix, uses it 
in determining ee CntmeticMeE ca he Wldcec 82° 25 G2ven if 
Appendix C. The program listing is also included. 
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The conventicnal beam pattern is determined by using 
equation z.£1 where rat is calculated by approximating the 
plot of relétive travel time vs. depth by a straight line at 
the receiving hydrophone depths. For ¢xample Figurs 3.3 
represents this plot for the 380 meter source depth. The 
Straight line is determined by a least squares linear 
regressicn which mirimizes the sum of the squares of the 
deviaticns of the actual data points from th? straight line 
Oeepest fit. Note that only data points which are on the 
dominant curve are used in calculating <he SEarghe lane. 
From the sgtraight line, relativ2 travel tines to the 
receiving hydrophones are calculated. The relative travel 
times for the 380 meter source depth are given in Table IL. 
They ccrrespond toa plane wave arrival angis of 3.73 
degrees. Ty) is determined by convertin these relative 
travel times +o phase delays using equation 2.12 and «hen 
taking the appropriate sine and cosine values of thes 
delays as in equation 2.24. The amplitude waights are 
initially assumed to be unity. 

A second method for obtaining the conventional bean 
pattern is calculated by the same procedure excep the 
amplitude weights which are determined by equation 2.50, the 


"9 
mL 
hydrorfhone. 


‘ amplizcude amplitude weights, are applied to each 
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hee a 


TABLE I 
Relative Travel Times For 380 Meter Source Depth 


Hydrophene Lepth Relative Tzavei Tine 
100 meters =U20) a055095 Seq. 
z00 meters -0.07131599 sec. 
200 meters -0.06697690 sac. 
400 meters -0.06263780 séc. 
500 meters —UmOge270? I sac. 
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A. EXACT SOLUTION 


An exact solution is derived for the beam pattern if the 
number of receiving hydrophones equal the numb2r of scurce 
depths. For examples, when the 5 receivers are used to 
discriminate between 5 source depths (220, 240, 260, 280, 
and 300 meters) there are 10 equations with 10 unknowns. 
Figure 4.1 is a plot cf the resulting bean Pate Goes ch ene 
source at the shallowest depth. Note that because of roun- 
doff errors in the *IMSL' subroutines there is a smali value 
for the resulting beam pattern at the non-energy source 


depths under investigation. 


B. FOUR DEPTHS WITH TWO RECEIVERS 


mem SOULcCS depths at 220, 280, 3090, and 320 meters with 
receiving hydrophones at 100 and 200 meters, there are 8 
equations with 4 unknowns. EXguce 4,2 is a plot c 
resulting Leam pattern with a source at the 220 meter depth. 
Figure 4.3 is the plot for the source at the 280 meter 
depth. 


C. TWENTY DEPTHS WITH TWO RECEIVERS 


For all 20 source depths with receiving hydrophones at 
100 and z00 meters, theres are 40 equations with 4 unknowns. 
Figure 4.4 is a plot of the resulting b2am pattern for 2 
souzce at the 220 meter depth. Figure 4.5 is a plet for the 
source at the 380 meter depth. 

In order to determine if the 1gt ealegia-cd in this case 


is the best an alternate method is devised. Four source 
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depths (220, 240, 260 meters, and another source test d 
with +¢he criginal source at the 220 meter depth and 
rece-vers at 100 and 200 maters are used. The beam pa 
is calculated each time with a different source test depth 
substituted for the fourth source depth. The 5 be 
resulting beam patterns are selected along with the 8 source 
depths (220, 240, 260 metézrs, and the 5 test depths wh 
created tke 5 best team patterns). Then, using these 8 
source deérths, igt is determined for the two receivers. 
This 1gt is applied to ths ‘two receivers and the bsan 
mencerp Cktained for all 20 source depths. 

The resulting beam pattern obtained by this alceznate 
methcd isn't as good as the beam pattern ob~ained by using 


eueez0esource depths in the determination of ‘'e'. 


D. TWENTY DEPTHS WITH FIVE RECEIVERS 


Fer all 20 source depths with all 5 receiving hydro- 
phones there are 40 equations with 10 unknowns. Figure 4.6 
is a plot cf the resulting beam pattern with ‘the source at 
the 220 meter depth. Figures 4.7, 4.8, and 4.9 are the 
plots for the source at the 360, 380, and 400 meczer depths 
respectively. 


E. CONVENTIONAL BEAMFORMER 


Figure 4.10 is a plot of the beam pattern for a convern- 
tional beamformer using linéar phase shifts across the array 
with the scurce at the 380 meter depth and the amplitude 
weights set to unity. All 20 source depths and 5 receiving 
hydrophones are used. Note that in figure 4.10 that *here 
is less than 1 db discrimation between each of the source 
depths. Figure 4.11 is the plot obtained for the amplitude 


weights set to values determined by equation 2.50. 
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F. RANGE OF 250 KILCMETERS 


The calculations were repeated for a range cf 250 kno 
using the same 20 Source depths and 5 receiving hydrophecnes. 
Figures 4.12 and 4.13 represent the plots cf relative <ravel 
times versus depth fcr the 220 and 380 meter source depths 
respectively. Figure 4.14 represents the straight line 
approximation of the relative travel times for the 380 meter 
depth. Note that in this figure the relative travel times 
are represented by twe straight lines; the upper line repre- 
sents linear phase shifts of the slowar travel times for che 
conventicnal beamformer while the lower lint tepresents 
linzar phase shifts of the faster travel times. Tables II 
and III are the straight line interpolations of zhese slower 
and faster travel +imes which correspond to arrival angles 
of 4.04 and -3.16 degrees respectively. Figures 4.15 and 
4.16 are the resulting beam patterns for the conventional 
beamfcrmer for the slewer travel times using unity amplitude 
weights and linear minimun variance amplitude weights 
respectively. Figures 4.17 and 4.18 are the beam patterns 
for the faster travel times. 

Figures 4.19, 4.20, 4.21, and 4.22 represent plots cf 
the beam pattern for the source at the 220, 340, 350, and 


380 meter depths respectively. 
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lS lS j 
| TABLE ITI 
Slower Ray Travel Times For 250 km Range and 380 a | 
| Source 
Hydrophone Lepth Relative Iravei Time 
|! 100 meters -0.09949109 sec. 
| z00 meters -0.09478615 ssc. | 
| 200 neters -0.09008%421 sec. 
400 meters -0.08537627 sec. | 
| 500 meters =). 08067 133 Ssec7 ! 
| | 


| 
| 
| 


4 
| 
| 
| 
| 

= 


TABLE IIL 
Faster Travel Times For 250 km Range and 380 m Scurce 


| | 
| 

| | 
| ' | 
| Hydrophene Depth Retaetve) Ieavele tne | 
100 neters -0.10100745 sec. 

| z00 meters -0.10468763 ssc. 
300 meters -0.10836781 sec. 
| 400 meters -0.11204799 ssc. 
| £00 meters -0.11572817 sec. | 
eee ee CO 
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V. DISCUSSION OF RESULTS AND RECOMMENDATIONS 

For many of the cases examined, the linear minimum vari- 
ance estimation technique gives high résoluticon in the 
“depth beam pattern" for sources at long range. Figures 4.2 
and 4.3 represent beam patterns for four depths with two 
receivers where the rain lobe is at the desired source depth 
with a beamwidth of 30 meters. Poo 620 depths with 2 
receivers the system is highly overdatermined and the beam 
pattern portrayed in figure 4.4 has its main lobe 80 meters 
from the scurce depth. At the source depth the strength of 
the beam pattern is 3.6 db down. PGE ti gure 4.5 th>= mein 
lobe is on the source depth with a beamwidth of 250 meters. 

For all 20 source depths with 5 hydzsophone rec¢sivers the 
results range from a beam pattern having its main lobe on 
the source depth (figure 4.7) with a beamwidth of 20 mnetears 
and a sécendary lobe 12.1 db down to a beam pattern having 
its main lebe 40 meters from the source depth (figure 4.22). 
For the Leam pattern in figure 4.22 the width of the main 
lobe is 80 meters and the strength of the beam pat=ern at 
the scurce depth is 0.8 db down. 


Ww 


=] 


When using the conventional beamformer, the beam patter 


rf 


a 


@ 


results with amplitude weights determined by the lin 
minimum variance estimation method were superior +o the bean 
pattern determined by using unity amplitud¢ weights. 
However in all cases the conventional beamfcrmer was sigrnif- 
icantly inferior to the beam pattern determined by the 
linear minimum variance estimation technique in terms of 
depth resolution. 

The linear minimum variance estimation technique could 
not be used for all ranges with the chosen sound speed 
eEotile. For ranges of 187 km and 235 km the relative 
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travel times preduce an ‘A' matrix Whee leeks sS SO i= 
conditioned that the ‘IMSL' subroutine ‘LINV2F' cannot 
determine an accurate inverse? of 'A 

There is sufficient evidence frem this initial investi- 
gaticn that the linear minimum variance estimation technique 
applied to a long linear vertical array can yield high reso- 
lution depth information about passiv= sources at very lcng 
ranges. However, further investigation is needed before any 
Peela scests are in order. Recommendations for further 
research are: 

1. The use of a more realistic speed of scund profile, 
Presctably One detuadly characteristic of the "SCFAR' 
channel. 

2. The use cf mere hydrophones in the vertical array. 
As more receivers are used the beam pattern should 
approach the desired beam pattern more closely (a 
less overdetermined system should yield smaller «octal 
tinimum méan squaré error). 

3. Investigation into causes for the peak response 
falling at other than the desired depth and techni- 
ques for correcting this problem when using the 
linear minimum variance estimation techinique. 

4. Alternate assumption fer choosing the ray paths <9 
ine iuds in the “A matrix. For exemple cne might 
choose the rays with the greatest intensity instead 
of the shortest travel time. 

5. Tke possibility of estimating range as well as th 
depth 9f a passive source with linear minimum vari- 
ance estimaticn techniques. 

In conclusion, the linear minimum variance estimation 
technique of beamforming was significantly superior to the 
conventional beamformer. HrgiesesoiN ton depth informacion 


about passive sources at long ranges is provided. 
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RELATIVE TRAVEL TIME CALCULATION 


This program calculates the relative travel times from 
Seem source depth tc the herizontal range and the depth of 
the ray at this range. To save paper, only the two shal- 


lowest scurce depths are listed. 
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INTERPOIATION OF RELATIVE TRAVEL TIME CALCULATIONS 
Qhis f~rograg does an interpolation of the ouput data 


generated in Appendix A. The relative travel times are 


interpolated for the receiving hydrophone depths. 
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APPENDIX C 
RESULTING BEAM PATTERN FOR CALCULATED WEIGHTS 


This program uses the output of Appendix B to calculate 
Boose A" Matrix. From this, the amplitude and phase weights 
are determined by the linear minumum variance? estimation 
technique. These weights are applied to the array and the 


beam pattern is obtained. 
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